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Abstract 



A factorized and analytical form for the A(7,pp) and A(e,e'pp) cross section 
is proposed. In the suggested scheme the two-proton knockout cross sections 
can be directly analyzed in terms of the ground-state correlation functions. 
Central, spin-spin and tensor correlations are considered. In the longitudinal 
channel, the (e,e'pp) cross section is predicted to exhibit a peculiar sensitivity 
to ground-state correlation effects. 
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The nuclear shell-model is a well-established theory for understanding the structure of 
atomic nuclei. Despite its obvious success, for a long time it has been realized that the strong 
repulsive nature of the nuclear force at short internucleon distances is likely to be at the 
origin of nuclear effects that are incompatible with the independent-particle nature of the 
shell-model. Obvious signs for short-range effects have recently been obtained from quasi- 
elastic (e,e'p) reactions. Extensive programs at several high-duty electron facilities predict 
values for the spectroscopic strength in the low-energy part of the residual-nucleus spectrum 
that are consistently lower than what was expected in the independent-particle model (IPM). 
This deviation has been the subject of several investigations, most of them pointing to short- 
and long-range correlations to be at the origin of the "missing" spectroscopic strength [|T],3 . 
In analyzing and interpreting the quasi-elastic (e,e'p) data it turned out to be extremely 
useful to make use of the following factorized form of the differential cross section || : 

d^<y 

de'dnJn p dTr EpaeppiWm ^ Exh (1) 

where E (p) is the energy (momentum) of the detected proton and a ep the elementary cross 
section for electron scattering on an off-shell proton. The spectral function p(\p m \ , E x ) 
is related to the probability of removing a nucleon with momentum p m from the target 
nucleus and finding the residual nucleus at an excitation energy E x . Strictly speaking, the 
above factorized form for the (e,e'p) cross section is only valid in the plane wave impulse 
approximation (PWIA) [[|, which puts aside effects like final state interactions (FSI) and 
photon absorption on two-body currents. 

A better understanding of the short-range correlations in nuclei is believed to come from 
a profound study of (7,pp) and (e,e'pp) reactions. The underlying idea is that reactions 
with an electromagnetic probe which induce two particles to escape are extremely sensitive 
to the two-body dynamics inside the target nucleus. When the main purpose of the two- 
nucleon knockout studies is the short-range part of the two-nucleon dynamics, proton-proton 
knockout is to be preferred above proton- neutron knockout. Indeed, the latter are also 
sensitive to correlations of the charge-exchange type (the most important one being one-pion 
exchange) that can heavily mask the effects of short-range nature Q. With two protons in 
the final state one has an enormous freedom when it comes to determining the position 
of the hadron detectors. Therefore, one could benefit from a simplified form for the cross 
section along the lines of Eq. ([!]). For this simplified form to be useful, one expects it to have 
some predictive power so that it can used to map the main sensitivities of the cross sections 
and optimize the kinematical conditions. In this paper we aim at deriving such a factorized 
form for the (e,e'pp) and (7,pp) cross section. Factorized expressions for the (7,pn) reaction 
have been derived in Refs. |5|||. The effect of factorization was studied in Ref. [[?]] and found 
to be a reasonable approximation at higher photon energies (u >150 MeV). Therefore, at 
moderate and larger values of the momentum transfer, it is to be expected that a factorized 
form for the (e,e'pp) cross section provides reasonable estimates. 

In order to arrive at a factorized form of a coincidence differential cross section it is a 
common procedure to work in a plane-wave model for the outgoing nucleons. Furthermore, 
we adopt the spectator approximation which means that the residual A-2 nucleons are 
assumed not to participate actively in the reaction process. Similar assumptions are at the 
basis of the factorized form ([I]) for the (e,e'p) cross section. In the two-nucleon knockout case, 
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however, the above assumptions are not sufficient for the cross section to factorize. This can 
best be explained by considering the two-nucleon knockout process as schematically depicted 
in Fig. |I[ There we assume that the photon couples to a correlated pair of protons. The 
center-of-mass (com) momentum P of the pair can be straightforwardly derived and reads 
P — ki t i+k2,i = ki + k2 — q- The relative momentum p re i = ( ki,i — k 2 ,ij /2 on the other hand, 
will generally depend on the momentum pc that is exchanged between the two correlated 
nucleons. Thus, in its most general form the cross section involves an integral over the pc 
which prevents it from being written in a factorized form. Now, the key approximation 
for the two-nucleon knockout cross section to factorize is that pc ~ 0. This assumption 
is related to the "quasi-deuteron approximation" which adopts the view that dinucleons in 
finite nuclei are mainly residing in relative S states ]7|||. Let pc = 0, then the relative 
momentum of the pair is determined by 

-* k\ — &2 q 

Prel = k± = ± - , (2) 

respectively corresponding with the situation that the photon is absorbed on proton "1" 
(k ) and "2" (k + ). With all of the above assumptions the (e,e'pp) cross section can be 
written in a factorized form : 

d^ 

E 1 p 1 E 2 p 2 crepp(k + ,k_,q) F hlM (P) , (3) 
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with k± =\ k± \ and F hlfl2 (P) denoting the probability to find a nucleon pair in the single- 
particle states hi(nili) and ^2(^2^2) with center-of-mass (com) momentum P. In the IPM 
this function reads : 

Fh 1 te(P)= E I dRe ip ^ nihm Amn 2 i 2m AR)\ (4) 

with (pnim the single-particle wave function for the orbital (nlm). The elementary cross 
section a epp can be considered as the equivalent of the a ep in Eq. (P and describes the 
physics of virtual photoabsorption on a diproton embedded in the target nucleus. Generally, 
a epp will depend on the photoabsorption mechanisms and the relative motion of the pair. It 
is precisely in the relative motion that the highest sensitivity to the correlation effects could 
be expected. We outline a method to derive an analytical expression for a epp . In doing this, 
we account for dinucleon correlations that go beyond the IPM. 

Adopting the PWA for the outgoing particles one has to evaluate the following type of 
matrix elements when calculating two-nucleon knockout cross sections 



M^ 2 m J « M « = Jdn...J df A A [e-^e-^ 2 ^ 
x ( l/2m Sl , l/2m S2 



h 2 jRM R 
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E4 1] ^+ E Jf{h3) \ Mh2,...,A)) , (5) 

i=l i<j=l 

where A is the anti-symmetrization operator and m Sl (m S2 ) the spin of escaping particle 
1 (2). The wave function for the residual nucleus ^^_2 Jr r is the two- hole state that is 
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created after knocking two protons out of the target nucleus. The operators jW and Jj$ 
are the one- and two-body parts of the nuclear current. As proton-proton emission is not 
sensitive to charge-exchange, the major component of the two-body operator is the isobaric 
current. The J^l is related to the charge density operator J2i Ge {q^Qh) e«5(r — fty. The 
transverse components of the one-body current JW (^ = ±1) is determined by the convection 
and magnetization current. The target wave function in the above expression is written as : 

^(1, 2, A) = F(l, 2, A)^(l, 2, A) , (6) 

where ipi is the IPM wave function and the operator F induces the correlations. In general, 
the operator F has many components, reflecting the full complexity of the nucleon-nucleon 
interaction. Calculations, however, have shown that the major correlation effects can be 
incorporated by considering an operator of the form PJIDj : 

A 

F(1,2,...A) = S J] [fc{r lj ) + Uar{r l3 )a l .a 3 + ftAr l3 )S lJ )f l .T J } . (7) 
i<j=i 

The first term accounts for central short-range correlations (commonly referred to as Jastrow 
correlations), whereas the other two induce spin-spin and tensor (SV,) correlations. 

The operator S is the symmetrization operator. 

In determining a epp we sum over the spins of the escaping nucleons. Further, instead 
of considering the contribution to each individual state we evaluate the cross section for 
emission out of a particular shell-model combination (niZ^n^a)- This means that we 
compute the averaged integrated cross section for a range of excitation energies in the A-2 
system. The range of excitation energies will be rather narrow when hole states close to the 
Fermi level are probed and grow wider as one or both nucleons are escaping from a deeper 
lying shell. 

One can expand the matrix element (|) in terms of g = 1 — fc, far and f tT . In doing 
this we have adopted the so-called "single-pair approximation" (SPA) ||11|| . This procedure 
is equivalent with an expansion into first order in g, f aT and f tT , retaining only those terms 
that contain the coordinates of both active nucleons. The SPA was earlier applied in Ref. 
in the context of (7,pn) reactions. Physically, the SPA is equivalent with multiplying the 
IPM relative wave function of the active pair with the correlation operator of Eq. [7|. With 
all these assumptions one obains the following general expression for a epp : 



epp ^Adfrec 



it f-ql 2@e\ ll 1 9«/ f\ Oe 

7± wl + yw +tan ^) WT+ 2^ WTT + 72V (e + e )tan ^ WLT 



(8) 



where f rec is the recoil factor and <jm the Mott cross section. The w's depend on the current 
operator and the different terms in the correlation operator (|7|). After lengthy calculations 
one arrives at the following analytical expressions : 

w L = 4e 2 (g(k+) + g{kJ)f {G E {q^)f + 40e 2 (f aT (k + ) + U T {k-)f {G E {q^)f 
+24e 2 (g(k + ) + g(k-)) (U T (k+) + /„(*_)) {G E (q^)f 

+ f /|e 2 (g(k + ) + g(k-)) (f? T (-k + ) + fl(-k-)) {G E {q^)f 
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It should be noted that the tensor and spin-spin correlations have only been considered for 
the longitudinal channel wl- We have defined the xz plane as the electron scattering plane 
with q along the z axis. In deriving the isobaric current contribution to the above expression 
we have used the operator specified in Ref. |T^] from which also all the coupling constants are 
taken. The G& is the A propagator in which an energy- dependent A 33 has been introduced 
I"2|| . The functions g(k), f aT {k) and ft T {k) occuring in the above expression are the Fourier 



transforms of the central functions occurring in the correlation operator (|7|) : 

g(p) = J dfe^(l - f c (r)) 
f aT (p) = J dfe^f aT (r) 
fl{p) ee / dfe^ r Y 20 (Q)f tT (r) . 



(10) 



In the absence of ground-state correlations only the transverse A current would give a non- 
vanishing contribution to a epp . The longitudinal channel Wl is totally determined by the 
ground-state correlations. 

Note that the a epp as written in Eq. (|8]) has formally the same form as the a ep in the 
(e,e'p) case. The physics situation is, however, very different. In the (e,e'p) case, the a ep 
contains information on electron scattering on a bound, off-shell nucleon. In quasi-elastic 
kinematics, this scattering process is predominantly sensitive to the one-body aspects of the 
target nucleus which are relatively well understood. For the (e,e'pp) case the situation is 
completely different. Here, a epp contains information on the two-body aspects of the nuclear 
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system. The Fh lt h 2 (P) occuring in the factorized (e,e'pp) cross section (^|) is a less- challenging 
quantity. As it deals with the c.o.m. motion of dinucleons it is rather insensitive to nucleon- 
nucleon correlations at small and moderate values of P. Recent calculations |nj predict 
some sensitivity at large momenta P. However, by then the function Fh ly h 2 (P) becomes so 
small in absolute magnitude that extremely small two-nucleon knockout cross sections can 
be expected. More favourable kinematical situations for probing correlations are created 
when considering relatively small values of P : in those cases larger cross sections are faced 
and the IPM predictions for F hljfl2 (P) can be considered as realistic choices. 
The (jiPp) cross section is uniquely determined by the wt term and reads : 



d 3 



a 



dQidQ2dT P2 

with, 



E 1 piE 2 p2cr 1P p (k+, k-, q-y) F hlM (P) , (11) 



a WP (k+, fc_, g 7 ) = 2E (2tt) 5 ^ Wt ' ^ 

In the transverse response function wt the contribution from the isobaric current and the 
terms related to ground-state correlations are competing. The Q 2 = —q^ dependence of 
the different terms in wt is investigated in Fig. |2|. We have considered in-plane kinematics 
(which means that both protons are escaping in the electron scattering plane) with a fixed P 
(which is chosen to point along the x axis), d\ and <pi (polar and azimuthal angle of escaping 
proton 1). For the results of Fig. ||] the central correlation function of the Gaussian type 
g(r) = 0.51e 1 " 52 ^ m > r as suggested in Ref. (13| was used. For all three energy transfers 



considered, the A (1232) current produces the largest contribution at the real photon point 
(Q 2 =0). Accordingly, (7,pp) reactions are predicted to exhibit a rather small sensitivity to 
Jastrow correlations. The second term in Wt, which corresponds with photoabsorption on 
the one-body convection current, has a marginal effect on the cross section. The term related 
to the magnetization current (first term in the expression for wt) has a q 2 dependence which 
makes it a dominant contribution at higher Q 2 . 

In Fig. |^ the Q 2 dependence of the longitudinal term (wl) in a epp is displayed for the 
same kinematical conditions as for Fig. ^| The first and second term in wl are originating 
from the central and spin-spin correlations respectively. The interference between the spin- 
spin and central correlations gives rise to the third term of w^. We remark that within the 
adopted approximations the contribution from the tensor correlations to wl is restricted to 
an interference term with the central correlations. For the wl we have investigated the role 
of the different terms in the correlation operator (0). Whereas central Jastrow correlations 
are frequently addressed in the literature, studies that include also spin-spin and tensor 
correlations are far less abundant. Using variational techniques the functions fc, far and 
ft T have been determined starting from a realistic nucleon-nucleon interaction ||[10|. We 
have used the correlation functions from the 16 O calculations of Ref. ||. It is striking that 
the Jastrow correlation function fc that comes out of these variational methods is very 
soft. For the sake of reference, we show in Fig. |3] also the result when including only the 
Gaussian Jastrow correlation function that was introduced earlier. With this choice of the 
central correlation function the longitudinal (e,e'pp) strength is predicted to be considerably 
larger than the strength produced by the variational central correlation function fc- The 
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variational method, however, is a more complete model in the sense that it accounts for 
the full complexity of the ground-state correlations. It is clear from Fig. |3] that spin-spin 
correlations are predicted to contribute substantially to the longitudinal (e,e'pp) strength. 
The effect of tensor correlations is rather marginal. This is not too surprising as they are 



generally considered to be a predominant proton-neutron correlation effect [14]. Note that 



the predictions for wl depend dramatically on the model assumptions with respect to the 
ground-state correlations. 

In Fig. H we compare the predictions of the suggested factorized cross section with 
recent 12 C(e,e'pp) data ]T5| and with unfactorized calculations. Two different regions 
in the excitation spectrum of the A-2 system were considered : E x < 13 MeV and 
23 MeV< E x < 48 MeV. The first region can be attributed to (lp) 2 knockout and the 
second to both (lp)(ls) and (Is) 2 knockout. For both of these regions we have calculated 
the angular cross sections with the variational and Gaussian correlation functions. With 
both choices one gets cross sections that are compatible with the data. 

For the Gaussian results we have compared the predictions of the factorized approach to 
the results of an unfactorized calculation. These calculations are based on an extension of 
the model outlined in Refs. [pl|[17]] and involve the same physics components as the unfac- 
torized approach in the sense that both central correlations and A 33 effects are considered. 
Unlike in a factorized model, one is no longer bound to a plane wave description for the 
outgoing particles in an unfactorized approach Jl6|,|l8|,[l9| . The right panels of Fig. f| show 
the results of the unfactorized calculations with both plane and distorted outgoing proton 
waves. Consequently the sole difference between the plane wave unfactorized calculations 
and the analytical approach discussed above is the treatment of the relative motion of the 
initial pair. It is noted that the unfactorized calculations produce angular cross sections that 
are wider than the analytical predictions. This is not too surprising given the fact that the 
relative motion of the initial diproton is now handled in its full complexity. After all, the 
factorized predictions could be considered as reasonable given that they can be performed 
in just a fraction of the computing time that the unfactorized calculations consume. Even- 
though the analytical expression a epp should not be considered as a fully-fledged alternative 
for the cumbersome unfactorized calculations, it could help in optimizing the kinematical 
conditions and determining the major trends and sensitivities of the cross section. 

Summarizing, we have derived a factorized form for the (e,e'pp) and (7,pp) cross section. 
This implies rather simple analytical expressions that should give a more transparent handle 
on the different physics components of photoinduced two-proton knockout reactions. Within 
the method developed here we have illustrated the sensitivity of the (e,e'pp) cross sections 
to ground-state correlation effects. 
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FIGURES 




FIG. 1. The emission of two-protons from a nucleus via (virtual) photon absorption. 
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FIG. 2. The Q 2 dependence of the transverse term wt in a epp for three values of the energy 
transfer. The kinematical condition is fixed through : P=50 (MeV) l x , #i=135 ° and 4>i=0 °. 
The dotted (solid) line shows the contribution from photoabsorption on the one-body convection 
(magnetization) current. The contribution from intermediate A33 production is shown with the 
dashed line. In these calculations a Gaussian central correlation function was used. See text for 
further details 
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FIG. 3. As in Fig. [| but now for the longitudinal term wl- The dashed curves only include 
Jastrow correlations, the dotted curves include Jastrow and spin-spin correlations, and finally the 
solid curves include Jastrow, spin-spin and tensor correlations. All these results have been obtained 



with the variational correlation functions of Ref. [10]. The dot-dashed line includes only Jastrow 
correlations, but now calculated with the Gaussian correlation function of Ref. |R 
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FIG. 4. The 12 C(e,e'pp) cross section for e = 475 MeV, w=212 MeV and q=270 MeV/c. One of 
the proton scattering angles was fixed at 27°. The dotted line is the calculated contribution from 
the ground-state correlations and the solid line is the prediction when including both ground-state 
correlations and intermediate A33 creation. Results for two types of correlation functions are shown. 
The data are from Ref. fi~5|| . The dot-dashed (dashed) line shows the result of an unfactorized 
calculation with plane (distorted) outgoing proton waves. 
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